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FOREWORD 


This  report  describes  a  part  of  the  ongoing  effort  to  make  the  latest  version  of 
the  Naval  Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD)  aeroprediction 
code  more  applicable  to  the  high  Mach  number  range  often  associated  with  modern 
weapon  systems.  The  methods  presented  enable  the  computation  of  three-dimen¬ 
sional  heat  transfer  rates  and  recovery  temperatures  as  part  of  the  overall  high  Mach 
number  solution. 

The  work  described  in  this  report  was  supported  primarily  through  the  Office  of 
Naval  Research  (Dave  Siegel,  Code  213)  and,  more  specifically,  the  Surface-Launched 
Weapons  Technology  Block  Program  managed  at  NSWCDD  by  Robin  Staton.  Partial 
funding  for  comparisons  of  the  new  predictive  techniques  to  other  methods  was 
provided  by  the  Strategic  Defense  Project  Office.  These  funds  are  managed  at 
NSWCDD  by  George  Long,  G205.  in  support  of  program  leadership  by  Pete  Stafford, 
K105.  Appreciation  is  expressed  to  these  individuals. 

This  report  has  been  reviewed  by  Thomas  J.  Rice,  Head,  Aeromechanics  Branch 
and  Danny  L.  Brunson,  Head,  Missile  Systems  Division. 
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ABSTRACT 


Methods  have  been  incorporated  into  the  Naval  Surface  Warfare  Center, 
Dahlgren  Division  aeroprediction  code  to  permit  the  computation  of  heat  transfer 
rates  and  recovery  temperatures  as  part  of  the  high  Mach  number  solution.  A  mass 
balance  technique  has  been  included  to  determine  the  correct  boundary  layer  edge 
entropy  to  use  along  the  surface  of  blunt  bodies.  Refinements  were  also  made  to  the 
pressure  prediction  routines  to  remove  discontinuities  and  improve  overall  results. 
Comparisons  of  results  from  the  new  methods  with  those  from  more  advanced 
engineering  codes,  with  other  techniques  of  similar  technical  level  of  detail,  and  with 
experimental  data  show  good  agreement.  These  new  capabilities  make  possible  the 
rapid  computation  of  three-dimensional  heat  transfer  information  for  a  wide  range  of 
geometric  configurations  and  flight  conditions. 
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1.0  BACKGROUND 


The  accelerating  pace  of  technological  advances  in  recent  years  has  given  rise  to 
a  new  generation  of  potential  military  threats  for  which  appropriate  counter¬ 
measures  must  be  developed.  In  the  case  of  surface-launched  tactical  missiles,  the 
new  targets  that  must  be  intercepted  are  faster,  more  maneuverable,  and  more 
difficult  to  detect  and  track.  Many  design  and  operational  factors  must  be  considered 
if  successful  engagements  are  to  be  possible  under  these  conditions.  One  of  the  more 
critical  concerns  lies  in  the  severe  aerodynamic  environment  to  which  the  missile  is 
exposed  as  a  result  of  the  high  velocities  needed  to  deal  with  sophisticated  targets 
that  may  be  approaching  quite  rapidly  and  may  only  be  detected  at  relatively  close 
range.  Depending  on  the  given  situation,  effective  defensive  measures  could  call  for  a 
missile  capable  of  upper  supersonic  or  even  hypersonic  flight.  This  regime  is 
currently  being  investigated  by  the  U.S.  Navy  with  designs  under  consideration, 
such  as  the  antitactical  ballistic  missile  (ATBM)  and  lightweight  exo-atmospheric 
projectile  (LEAP),  which  call  for  upper  flight  Mach  numbers  of  around  15. 

As  the  performance  envelope  is  extended  into  the  hypersonic  region  above 
Mach  6,  two  important  aerodynamic  phenomena  become  increasingly  important. 
The  first  of  these  relates  to  the  manner  in  which  the  air  molecules  in  the  vicinity  of 
the  missile  store  energy  at  the  elevated  temperatures  created  by  the  strong  bow 
shock.  At  temperatures  up  to  about  1400  °R,  the  heat  dissipated  by  the  shock  is 
stored  in  the  translational  and  rotational  energy  modes  of  the  gas  molecules.  The 
fluid  behavior  is  adequately  described  by  the  perfect  gas  law,  and  the  specific  heats 
remain  constant.  As  the  temperature  rises,  the  vibrational  energy  modes  of  the 
molecules  become  excited.  ri  he  perfect  gas  law  is  no  longer  valid,  and  the  specific 
heats  become  functions  of  temperature.  Above  about  4500  °R,  the  air  molecules  begin 
to  dissociate  and  chemical  reactions  occur  leading  to  even  greater  complexities.  This 
departure  from  perfect  gas  behavior  is  often  referred  to  generically  as  real-gas  effects. 
It  tends  to  become  significant  above  approximately  Mach  6.  The  overall  result  is  to 
lower  the  gas  temperatures  below  those  that  would  be  predicted  by  a  perfect  gas 
analysis.  The  increase  in  modes  of  internal  energy  storage  mean  that  less  energy 
goes  into  the  translational  mode,  which  is  directly  measured  as  temperature.  This 
effect  will  not  be  discussed  in  detail  in  this  report  other  than  in  the  context  of  the  use 
of  appropriate  real-gas  relations  in  high-temperature  computations.  For  further 
information,  Reference  1  provides  a  general  discussion  and  Reference  2  explains  in 
detail  how  the  aeroprediction  code  has  been  modified  to  model  real-gas  physics. 


1-1 


NSWCDD/TR-93/29 


The  phenomenon  described  is  characteristic  of  all  high  Mach  number  flows  and 
is  unrelated  to  the  fluid  viscosity.  Thus,  real-gas  effects  must  be  included  even  in 
inviscid  flow  models  such  as  the  Euler  equations.  The  second  effect  that  becomes 
important  in  high  Mach  number  flows  is  related  directly  to  the  viscosity  of  the  fluid. 
Boundary  layer  theory  tells  us  that  the  fluid  particles  in  contact  with  the  body 
surface  are  at  rest  relative  to  it.  This  is  the  viscous  no-slip  condition.  However,  in 
the  case  of  supersonic  or  hypersonic  flight,  the  particles  at  some  distance  above  the 
surface  are  moving  at  high  velocities.  Thus,  there  is  a  region  separating  the  body 
surface  and  the  external  flow  within  which  the  fluid  particles  are  decelerated  from 
some  large  velocity  to  rest.  This  is  the  classical  boundary  layer,  and  it  is  usually 
quite  thin,  although  in  some  hypersonic  flows,  it  can  become  more  pronounced.  The 
boundary  layer  is  dominated  by  viscous  shear  forces  that  dissipate  the  kinetic  energy 
of  fluid  particles  approaching  the  surface.  Energy  cannot  be  destroyed,  so  this 
process  results  in  an  increase  in  the  internal  energy  of  the  fluid,  the  magnitude  of 
which  is  directly  related  to  the  kinetic  energy  of  the  external  flow.  This  increase  in 
internal  energy  appears  in  the  form  of  a  rise  in  boundary  layer  temperature  and, 
because  kinetic  energy  is  proportional  to  the  square  of  the  velocity,  it  is  apparent  that 
as  the  Mach  number  increases,  the  boundary  layer  temperature  can  increase  rapidly, 
leading  to  greater  transfer  of  heat  to  the  body  surface.  The  Mach  number  at  which 
this  boundary  layer  heating  becomes  an  important  factor  is  dependent  on  a  number  of 
factors,  such  as  the  composition  and  structure  of  the  missile,  time  of  flight,  internal 
environmental  constraints,  etc.  In  some  situations,  it  can  become  significant  at 
supersonic  velocities.  At  hypersonic  Mach  numbers,  it  cannot  be  ignored  under  any 
circumstances. 

This  report  describes  the  methodology  that  has  been  incorporated  into  the 
aeroprediction  code  to  provide  estimates  of  the  boundary  layer  heating.  Comparisons 
are  presented  with  results  from  an  approximate,  point-solution  aeroheating  code 
(MINIVER)  and  several  more  sophisticated  viscous  flow  solvers,  as  well  as  with 
experimental  results.  A  comparison  with  inviscid  flow  properties  computed  by 
MINIVER  is  also  presented  to  validate  the  real-gas  techniques  previously 
incorporated  into  the  aeroprediction  code. 
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2.0  INTRODUCTION 


The  task  of  generating  numerical  predictions  of  boundary  layer  heating  can  be 
approach  ;d  in  several  ways.  The  most  rigorous  treatment  involves  a  solution  of  the 
complete  time-dependent  Navier-Stokes  equations  including  chemically  reacting 
species  if  the  temperatures  are  sufficiently  high.  Given  present  computer  resources, 
this  approach  is  much  too  costly  to  be  used  as  a  general  engineering  tool,  requiring 
several  hours  of  computational  time  and  several  megawords  of  memory  on  a  super¬ 
computer.  It  would  have  to  be  reserved  for  a  select  number  of  special  cases  with 
complicating  factors,  such  as  flow  separations  or  shock-boundary  layer  interactions, 
which  are  not  amenable  to  simpler  solution  tactics. 

If  the  full  equations  are  simplified  by  dropping  the  streamwise  (and,  in  some 
formulations,  the  crossflow)  viscous  terms,  the  resulting  Thin  Layer  Navier-Stokes 
equations  are  somewhat  easier  to  handle  while  still  giving  excellent  results  for  most 
external  flows.  However,  the  solutions  are  costly,  requiring  only  slightly  less 
computer  time  and  storage  than  the  complete  equations  and  so  are  still  not  suitable 
for  routine  engineering  use. 

Further  simplification  can  be  achieved  by  dropping  the  time-dependent  terms  to 
produce  the  Parabolized  Navier-Stokes  (PNS)  equations.  Obviously,  these  are  valid 
only  for  steady  flows,  but  they  achieve  great  savings  in  computer  resources  by 
allowing  space-marched  solutions  that  require  storage  of  only  two  or  three  crossflow 
planes  of  the  computational  grid  at  once  instead  of  the  entire  field.  This  method 
cannot  be  applied  if  there  are  streamwise  separations  or  regions  where  the  flow 
exterior  to  the  boundary  sr  is  subsonic.  The  PNS  equations  can  be  a  useful  tool 
during  the  design  process,  but  with  computation  times  measured  in  minutes  or  even 
hours,  the  cost  still  precludes  their  use  in  the  preliminary  phase  to  generate  large 
matrices  of  data  for  multiple  configurations. 

With  a  final  simplification  of  dropping  the  viscous  terms  entirely,  the  widely 
used  Euler  equations  are  obtained.  If  the  time-dependent  terms  are  retained,  it  is 
possible  to  model  unsteady  and  subsonic  flows.  Without  these  terms,  only  steady 
supersonic  flows  can  be  handled,  but  it  becomes  possible  to  use  a  rapid  space-marched 
solution  technique.  If  the  interest  is  only  in  the  surface  pressures  on  the  given 
missile,  these  equations  can  provide  good  solutions,  but  if  viscous  heating  effects  are 
important,  the  inviscid  Euler  results  must  be  coupled  to  a  boundary  layer  solver  of 
some  type.  This  approach  is  possible  because  of  the  relative  weakness  of  the  viscous- 
inviscid  interactions  in  many  problems.  To  a  first  approximation,  there  is  no 
pressure  change  across  the  boundary  layer  allowing  the  use  of  the  Euler  surface 
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pressures  for  all  boundary  layer  computations.  In  addition,  at  sufficiently  high 
Reynolds  numbers,  the  boundary  layer  thickness  is  small  in  relation  to  the  dimen¬ 
sions  of  the  missile  body.  Thus,  the  viscous  surface  layer  has  no  significant  effect  on 
the  shape  of  the  body  as  seen  by  the  outer  inviscid  flow.  The  Euler  equations  can  be 
°olved  to  get  the  outer  flow  field  and  the  computed  surface  pressures  and  velocities 
can  be  used  as  upper  edge  conditions  for  the  boundary  layer  solution.  Even  if  the 
boundary  layer  thickness  increases  to  the  point  that  it  may  no  longer  be  considered 
thin,  it  may  still  be  possible  to  use  this  coupled  approach  by  employing  an  iterative 
sequence  in  which 

1.  Inviscid  flow  is  computed  with  the  Euler  equations. 

2.  Boundary  layer  equations  are  solved  using  the  inviscid  surface 
conditions  as  inputs  giving,  among  other  things,  the  boundary 
layer  thickness. 

3.  Adjustments  are  made  to  the  body  shape  reflecting  the  displace¬ 
ment  effect  of  the  boundary  layer. 

4.  The  sequence  is  repeated  until  convergence  is  achieved. 

At  higher  hypersonic  Mach  numbers,  the  boundary  layer  can  occupy  a  significant 
portion  of  the  space  between  the  body  and  the  shock  or,  in  some  cases,  the  outer  edge 
of  the  boundary  layer  may  merge  with  the  shock.  When  the  viscous-inviscid  inter¬ 
action  becomes  this  strong,  even  the  iterative  approach  no  longer  gives  good  results, 
and  one  of  the  previously  discussed  methods  would  have  to  be  used. 

The  Euler  equations  are  used  on  a  fairly  routine  basis  at  all  stages  of  the  design 
and  engineering  process  to  generate  force  and  moment  data  for  aerodynamic  bodies. 
The  coupled  Euler/boundary  layer  approach  can  provide  good  estimates  of  surface 
heating  over  a  wide  range  of  conditions,  but  depending  on  the  particular  combination 
of  solvers  used,  it  can  require  computer  resources  on  the  order  of  those  used  for  the 
PNS  equations.  Using  space-marched  Euler  and  integral  boundary  layer  methods, 
computational  costs  can  be  reduced  to  the  point  that  engineering  applications  become 
possible,  but  a  typical  run  to  get  the  surface  heating  over  a  generic  missile  shape  for  a 
single  Mach  number-altitude  combination  can  still  require  computational  times  of 
several  minutes.  Thus,  this  approach  may  also  prove  to  be  unsuitable  for  use  in 
preliminary  design  investigation. 

Before  the  advent  of  modern  high-speed  supercomputers,  many  approximate 
analytical  techniques  were  developed  to  model  the  characteristics  of  a  wide  range  of 
specialized  aerodynamic  flows.  By  making  certain  simplifying  assumptions,  the 
complexity  of  the  governing  equations  was  reduced  to  a  level  that  could  be  handled 
easily  by  the  available  tools — often  nothing  more  than  hand  computations  being 
required.  This  report  will  not  attempt  to  present  a  survey  of  such  techniques  but  will 
make  reference  only  to  those  that  bear  directly  on  the  aeroprediction  viscous  heating 
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modifications.  It  is  important  to  point  out  that  such  techniques  exist  and  ti'at  they 
produce  acceptable  results  as  long  as  their  range  of  application  does  not  go  beyond  the 
restricting  conditions  inherent  in  their  development.  A  number  of  computer  codes 
have  been  put  together  using  such  techniques,  combining  several  analytical  methods 
and,  perhaps,  some  empirical  databases  to  cover  a  range  of  aerodynamic  conditions 
and  geometric  configurations.  The  major  advantage  of  this  approach  lies  in  its  speed, 
permitting  computations  in  seconds  or  even  fractions  of  a  second  on  any  mainframe 
computer  or  workstation.  Their  obvious  drawback  lies  in  their  limited  range  of 
validity,  but  these  limitations  often  are  not  severe,  and  a  large  percentage  of 
problems  likely  to  be  encountered  lend  themselves  to  such  treatment. 

The  aeroprediction  code3,4,5  is  a  collection  of  approximate  techniques  put 
together  by  Naval  Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD)  personnel 
to  generate  force  and  moment  data  for  a  wide  variety  of  geometric  configurations  over 
a  velocity  range  from  subsonic  to  Mach  20.  Limited  capabilities  are  also  available  for 
computing  skin  friction  drag,  but  no  viscous  heating  methods  have  previously  been 
incorporated.  Over  the  years,  it  has  proven  to  be  an  effective  design  tool  and  has  been 
distributed  to  approximately  50  private,  governmental,  and  international  clients. 6  It 
has  been  the  subject  of  several  upgrades  to  extend  its  range  of  application,  the  most 
recent  of  these  being  the  incorporation  of  high-temperature,  real-gas  effects  into  its 
second-order  shock  methodology,  thus  extending  its  usefulness  into  the  hypersonic 
velocity  regime. 2  As  pointed  out  previously,  at  these  high  Mach  numbers,  boundary 
layer  heating  becomes  a  dominant  factor  in  airframe  design  considerations.  Because 
this  was  not  an  area  that  had  been  addressed  in  previous  aeroprediction  code 
evolution,  it  was  decidjd  that  addition  of  this  capability  would  make  the  high  Mach 
number  version  a  much  more  useful  tool.  In  keeping  with  other  aeroprediction 
methodologies,  the  viscous  heating  computations  were  to  be  implemented  using 
approximate  analytical  techniques  requiring  a  minimum  of  computer  resources  to 
execute.  The  following  sections  describe  in  detail  how  this  goal  was  achieved. 
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3.0  HEAT  TRANSFER  COEFFICIENTS 


The  boundary  layer  heating  information  is  presented  in  the  aeroprediction 
output  in  the  form  of  a  heat  transfer  coefficient  and  an  adiabatic  wall  temperature 
(also  known  as  the  recovery  temperature)  at  each  computational  point.  These 
variables  are  related  in  the  following  manner. 


H  = 


(T  -T  ) 

aw  w 


(3-1) 


where  H  is  the  heat  transfer  coefficient,  qw  is  the  surface  heat  transfer  rate,  Taw  is 
the  adiabatic  wall  temperature,  and  Tw  is  the  wall  temperature.  The  adiabatic  wall 
temperature  is  defined  as 


T  =  T  +  r  (Ta  -  T  )  (3-2) 

aw  e  c  0  e 

where  Te  is  the  boundary  layer  edge  temperature  (inviscid  solution  surface 
temperature)  and  To  is  the  stagnation  temperature.  The  recovery  factor,  rc,  is 
defined  as 


r 

c 


(3-3) 


and  is  a  measure  of  the  degree  of  efficiency  involved  in  converting  the  kinetic  energy 
of  the  external  flow  into  heat  within  the  boundary  layer.  If  this  conversion  were 
100-percent  efficient,  the  adiabatic  wall  temperature  and  the  stagnation  temperature 
would  be  the  same.  Taw  is  the  highest  possible  temperature  that  can  be  reached 
within  the  boundary  layer  and,  by  extension,  this  is  the  highest  temperature  to  which 
the  wall  can  be  heated.  The  length  of  time  required  for  the  wall  temperature  to  reach 
this  upper  limit  is  dictated  by  the  rate  of  heat  transfer,  qw.  If  the  surface  actually 
reaches  Taw,  it  can  be  seen  from  Equation  3-1  that  the  heat  transfer  rate  would  go  to 
zero,  thus  the  term  adiabatic  wall  temperature. 

The  heat  transfer  coefficient  for  high  temperature  flows  is  often  expressed  in 
terms  of  enthalpies  as 


n  = 


(3-4) 
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where  x\  is  the  enthalpy -based  heat  transfer  coefficient,  haw  is  the  adiabatic  wall 
enthalpy,  and  hw  is  the  enthalpy  of  the  fluid  at  the  body  surface.  In  this  case,  haw  is 
defined  as 


=  h  +  r  (h.  -  h  ) 
e  c  0  e 


and  the  recovery  factor,  rc,  is 


r 

C 


(3-5) 


(3-6) 


At  temperatures  above  about  1500  °R,  the  enthalpy  formulation  is  the  more 
rigorously  correct  of  the  two.  Below  that  point,  they  are  equivalent,  because  the 
specific  heats  of  the  gas  are  constant  and  h  =  cpT.  Within  the  code,  all  real-gas,  high- 
temperature  computations  are  performed  using  the  enthalpy  as  the  basic  variable, 
but  to  maintain  consistency  with  other  codes  in  use  at  NSWCDD,  the  results  are 
converted  to  the  temperature  form  for  output.  In  either  case,  the  recovery  factor  can 
be  related  to  the  Prandtl  number  of  the  fluid  by 

r  =  V  Pr  (laminar)  (3-7) 

C 


and 


r  =(Pr)1/3  (turbulent) 

C 

The  heat  transfer  is  normalized  in  the  form  given  by  Equation  3-1  because  the 
coefficient  H  remains  fairly  constant  over  a  wide  range  of  wall  temperatures.  The 
heat  transfer  itself  varies  from  some  maximum  value  at  Tw  =  0  (cold  wall)  to  zero  at 
Tw  =  Taw  (adiabatic  wall).  Given  H  and  Taw  for  a  designated  Mach  number,  angle  of 
attack,  and  altitude  combination,  it  is  possible  to  compute  the  actual  heat  transfer 
rate  for  any  wall  temperature.  Then,  given  the  composition  and  structure  of  the  wall, 
its  rate  of  temperature  change  can  be  determined.  By  performing  these  computations 
sequentially  at  a  selected  number  of  points  along  a  missile’s  trajectory,  a  time  history 
of  the  surface  temperature  over  the  body  can  be  constructed. 
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4.0  BOUNDARY  LAYER  HEATING  METHODS 


The  user  has  the  option,  specified  in  the  input  file,  of  using  either  perfect-gas  or 
real-gas  relations  for  any  aeroprediction  run.  For  a  detailed  discussion  of  the 
theoretical  background  of  the  two  approaches  and  of  the  real-gas  extensions  applied 
to  the  second-order,  shock-expansion  theory  in  aeroprediction,  refer  to  Reference  2. 
The  important  point  to  bear  in  mind  is  that  as  flight  Mach  numbers  rise  beyond  about 
6,  real-gas  effects  become  increasingly  important.  Pressures  are  affected  to  a 
minimal  degree,  but  temperatures  are  lowered  substantially  at  the  higher  Mach 
numbers  as  new  modes  of  energy  storage  are  excited  in  the  gas.  In  the  discussion  that 
follows,  both  perfect-gas  and  real-gas  methods  will  be  set  forth.  Both  are  available 
for  use  in  the  code.  The  perfect-gas  option  gives 
applied,  but  it  is  not  recommended  above  Mach  6. 


4.1  NOSE  TIP  STAGNATION  POINT 

The  free-stream  flow  is  first  passed  through 
priate  aeroprediction  subroutines.  This  process 
methods,  depending  on  the  choice  of  input  options, 
fully  described  by  the  subroutine  outputs  that  are 
conditions  behind  the  shock) 

v2  =  Velocity 
p2  =  Density 
h2  =  Enthalpy 
e2  =  Internal  energy 
S2  =  Entropy 
a2  s  Sonic  velocity 
Y2  =  Specific  heat  ratio 
z2  -  Compressibility  factor 

Stagnation  conditions  behind  the  shock  are  needed,  and  these  can  be  obtained 
by  an  isentropic  compression  of  the  flow  to  zero  velocity.  The  stagnation  pressure  is 
given  by 


faster  run  times  when  it  can  be 


a  normal  shock  using  the  appro- 
uses  either  real-  or  perfect-gas 
The  flow  behind  the  shock  is  then 
(where  the  subscript  two  denotes 


Po  =  P2  + 


(4-1) 


4-1 


NSWCDD/TR-93/29 


This  relation  can  be  used  because  the  flow  behind  the  shock  is  subsonic.  In  fact, 
the  higher  the  free-stream  Mach  number,  the  more  subsonic  the  downstream  flow. 
Thus,  even  though  the  Bernoulli  relation  is  strictly  valid  only  for  incompressible 
flow,  the  low  Mach  numbers  behind  the  shock  minimize  any  compressibility  effects  as 
the  gas  is  decelerated. 7  Thus,  because  the  stagnation  pressure  can  be  computed  and 
the  stagnation  entropy  is  known  because  the  compression  is  isentropic,  the  remaining 
gas  properties  can  be  determined. 


4.1.1  Perfect-Gas  Stagnation  Point  Conditions 

The  specific  heat  ratio  is  assumed  to  be  constant  and  equal  to  1 .4.  Then 

V2"1 


T0  =  T2(i  + 


(4-2) 


Po  = 


c 

p 


R 


ho  cpTo 


M0  =  0 
zo  =  1 

R  is  the  gas  constant. 

4.1.2  Real-Gas  Stagnation  Point  Conditions 

For  real-gas  conditions,  first  compute  the  stagnation  enthalpy. 


(4-3) 

(4-4) 

(4-5) 

(4-6) 

(4-7) 

(4-8) 
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Then,  given  p0,  a  series  of  curve  fits  is  used  to  obtain  the  remaining  prop- 
erties.8,9 


P0  =  Po(po>ho) 


(4-10) 


T  =T. 


o(Po-Po) 


(4-11) 


eo  =  eo(Po*Po) 


(4-12) 


%  *Veo,po* 


(4-13) 


ao  =  ao(eo*so> 


(4-14) 


=  _Po_ 
Z°  PoRTo 


(4-15) 


Y°  RT„ 


(4-16) 


vo=0 


(4-17) 


M0  =  ° 


(4-18) 


4.1.3  Stagnation  Point  Heating 

Given  either  the  perfect-  or  real-gas  properties  established  in  the  preceding 
sections,  the  nose-tip  stagnation  point  heating  rate  can  be  found  froml 


qw=  0.763  Pr  06  V  -2  (h>w  -  hj 


(4-19) 
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Pr  is  the  Prandtl  number  and  is  considered  to  be  a  constant  0.7  for  the  perfect 
gas  case  and  is  found  from  the  curve  fits  for  a  real  gas. 


Pr  =  Pr  (Tff  p,) 


(4-20) 


The  coefficient  of  viscosity,  po>  is  determined  from  Sutherland’s  Law  for  perfect- 
gas  runs  but  is  computed  from  the  curve  fit  relations  for  real  gases. 


Pq  ^0  (To*  Po> 


(4-21) 


The  stagnation  point  streamwise  velocity  gradient,  due/dx,  can  be  expressed  as 

^  =  L  (4-22) 

dx  rN  Po 

where  a  spherical  nose  tip  of  radius  tn  is  assumed.  If  the  aeroprediction  input  file 
specifies  a  sharp  nose,  Tn  is  set  to  0.01  in.  The  adiabatic  wall  enthalpy,  haw,  is  defined 
in  Equation  3-5,  and  hw  is  the  enthalpy  of  the  fluid  at  the  wall  temperature.  There  is 
only  one  form  of  Equation  4-19  for  both  laminar  and  turbulent  flow.  At  the 
stagnation  point,  the  flow  will  always  be  laminar. 


4.2  CONTROL  SURFACE  LEADING  EDGE  STAGNATION  LINE 

It  is  assumed  that  control  surfaces  are  located  in  a  “  +  ”  configuration,  and  the 
boundary  layer  heating  is  computed  at  the  midspan  of  the  upper  and  lower  horizontal 
panels.  This  arrangement  is  consistent  with  that  presently  employed  in  the  high 
Mach  number  control  surface  subroutines  in  aeroprediction.  These  subroutines  use  a 
local  method  that  starts  from  free-stream  conditions,  and  then,  using  oblique  shocks 
or  Prandtl-Meyer  expansions,  turns  the  flow  through  the  necessary  angle  to  make  it 
parallel  to  the  surface.  No  interference  effects  are  included.  It  is  realized  that  this 
approach  can  be  improved,  and  future  work  will  be  devoted  to  this  area  as  circum¬ 
stances  permit. 

The  first  step  in  the  computations  involves  decomposing  the  free  stream  into 
components  normal  and  parallel  to  the  leading  edge. 

v  =  v  cos  A  (4-23) 

n  ® 


a/"2 - 2  (4-24) 

v  =  V  v  —  v 
p  »  n 
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The  n  and  p  subscripts  denote  directions  normal  and  parallel  to  the  leading 
edge,  respectively.  A  is  the  sweep  angle  with  A=0  corresponding  to  no  sweep.  The 
normal  velocity  component  is  passed  through  a  normal  shock,  which  produces  a  set  of 
downstream  conditions  as  given. 

vds.ii  =  Downstream  component  normal  to  leading  edge 

Tds  =  Downstream  temperature 

pos  =  Downstream  pressure 

pos  =  Downstream  density 

hos  =  Downstream  enthalpy 

eos  =  Downstream  internal  energy 

Sds.  =  Downstream  entropy 

aos  =  Downstream  sonic  velocity 

Yds  =  Downstream  specific  heat  ratio 

zqs  =  Downstream  compressibility  factor 

vp  is  not  affected  by  the  shock.  Assuming  that  the  leading  edge  is  analagous  to  an 
infinite  span  cylinder,  the  shock  wave  will  be  parallel  to  the  stagnation  line.  Because 
vp  is  parallel  to  the  shock,  it  will  also  be  parallel  to  the  stagnation  line  and,  similarly, 
VDs.n  will  be  normal  to  the  stagnation  line  and  must  go  to  zero  at  the  surface.  This 
loss  of  velocity  will  be  felt  as  an  increase  in  the  pressure.  Thus,  at  the  surface  the 
result  is 

v  '  =  v  (4-25) 


Vne  ->  0  (4-26) 

DS,n 


1  2 

PDS  -  PDS  +  2  PDSVDS,n  (4-27) 

Equation  4-27  is  a  valid  approximation  because  of  the  relatively  low  velocities 
involved,  as  argued  previously.  The  prime  superscript  denotes  conditions  at  the 
surface  where  the  normal  velocity  has  been  reduced  to  zero.  The  other  fluid 
properties  can  now  be  found. 


Again,  assume  a  constant  specific  heat  ratio  of  1.4,  then  proceed  as  follows: 


(4-28) 
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M  = 


DS,n 


DS 


YDS  1 


TDS  =  TDS (1  +  2  "D 


PDS 


DS 


RT. 


DS 


h  =  C  T 
DS  pDS 


ZDS  1 


aDS  =  v/  YDS  R  T  DS 


mds  = 


'DS 


DS 


,  -  PDS 

eDS  ”  hDS 

PDS 


SDS  "  SDS 


(4-29) 

(4-30) 

(4-31) 

(4-32) 

(4-33) 

(4-34) 

(4-35) 

(4-36) 

(4-37) 


In  addition  to  these  variables,  the  swept  leading  edge  heating  computations 
require  as  input  the  stagnation  conditions  behind  a  normal  shock  at  the  free-stream 
Mach  number.  These  are  obtained  in  the  same  manner  as  described  in  Section  4.1.1. 
It  is  also  necessary  to  determine  the  stagnation  conditions  behind  the  swept  leading 
edge  shock.  For  a  perfect  gas,  they  can  be  determined  as  follows: 
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P0  PDS  ^  + 


*DS 

^  «MV“- 


(4-38) 


<  •  ^DS  ^  >o 

T0  =  tds(1  +  T— “m* 


(4-39) 


Pf>  = 


RT„ 


ho  =cPTo 


V0  =  ° 


M0  =  0 


zo  =  1 


(4-40) 


(4-41) 


(4-42) 


(4-43) 


(4-44) 


4.2.2  Real-Gas  Stagnation  Line  Conditions 

The  real-gas  computations  begin  with  the  determination  of  the  enthalpy. 

hD8  =  hDS  +  2  VDSji 


(4-45) 


The  entropy  is  known  since  the  compression  to  the  prime  conditions  is 
isentropic.  Using  the  curve  fit  relations  once  again  gives 


PDS  “  PDS  (pD6’hDS^ 


TDS  "  TDS(PDS’  PD8) 


(4-46) 


(4-47) 
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aDS  aDS 'PDS’ 

(4-48) 

aDS 

(4-49) 

eDS  =  eDS  (pDS’  PDS^ 

(4-50) 

aDS 

VDS  rtds 

(4-51) 

ZDS  =  1 

(4-52) 

The  required  stagnation  conditions  behind  the  free-stream  normal  shock  are 
determined  as  in  Section  4.1.2.  It  is  not  as  straightforward  to  get  the  stagnation 
conditions  behind  the  leading  edge  oblique  shock.  The  Bernoulli  relation  cannot  be 
used  because  the  Mach  number  is  too  high  to  ignore  compressibility  effects.  The 
supersonic  relations  for  stagnation  conditions  cannot  be  used  because  they  do  not 
include  real-gas  effects.  An  isentropic  compression  to  the  known  stagnation  enthalpy 
using  real-gas  relations  is  required.  Unfortunately,  the  available  curve  fits  do  not 
contain  a  function  with  this  combination  of  variables.  It  was  necessary  to  employ  an 
iterative  approach.  First,  a  stagnation  pressure  was  assumed.  The  value  predicted 
by  the  supersonic  flow  formula  was  used  as  this  initial  guess.  Then,  using  the 
available  curve  fits,  a  corresponding  value  of  stagnation  entropy  was  computed  as 
follows: 

P0  P0  (p0> 

(4-53) 

eo  =  eo  Po) 

(4-54) 

S0  =S0(e0’P0) 

(4-55) 

The  entropy  value  was  checked  for  agreement  with  the  isen tropic  assumption. 
If  it  was  not  correct,  a  second  stagnation  pressure  value  was  chosen  and  the  process 
repeated.  If  the  resulting  entropy  was  still  off,  a  variable  secant  iteration  scheme  was 
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used  to  get  the  correct  stagnation  pressure.  The  remaining  flow  properties  could  then 
be  determined  using  the  curve  fit  relations. 


T0  T0  (p0,  p0) 


(4-56) 


ao  =  ao(Po-So) 


(4-57) 


(4-58) 


(4-59) 


v0=0 


(4-60) 


M0=0 


(4-61) 


4.2.3  Beckwith  and  Gallagher  Swept  Cylinder  Heating 

The  method  of  Beckwith  and  GallagherlO.il  is  used  to  compute  the  leading  edge 
heat  transfer  rate.  For  laminar  flow,  the  stagnation  line  heating  rate  is  given  by 


qw,i  =  0.57Pr~06 V  —  (h^  -  hj  (coaA)U 


(4-62) 


Equation  4-62  is  identical  to  Equation  4-19  with  the  exception  of  the  leading 
coefficient  and  the  addition  of  the  last  term  involving  the  cosine  of  the  sweep  angle. 
For  turbulent  flow,  the  relation  is  more  complicated. 


(a*u*)0,8  <  au  \0.2 

q  =1.04Pr-0'6 — — —  (vre,sinA)0'6(  — ^ )  (h  -h  ) 
Hw.t  ')0.6  DS  (  dx  /SL  «*  w 


du. \o^ 


(4-63) 


The  *  superscript  denotes  evaluation  at  the  Eckert  reference  enthalpyl2  which 
is  given  by 


h  =0.5(h  +  h  )  +  0.22  (h  -h  ) 

we  aw  e 


(4-64) 
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where  he  is  the  enthalpy  at  the  outer  edge  of  the  boundary  layer.  For  the  perfect-gas 
case,  it  is  acceptable  to  compute  a  reference  temperature  instead. 

T*  =  0.5  (T  +  T  )  +  0.22  (T  -  T  )  (4-65) 

we  aw  e 


The  reference  enthalpy  is  then  found  by  multiplying  T*  by  the  specific  heat,  cp.  The 
streamwise  velocity  gradient  at  the  stagnation  line  is  given  by 


du 

'dT'-. 


(4-66) 


where  r^  is  the  leading  edge  radius  of  the  control  surface  (which  is  assumed  to  be 
cylindrical).  If  a  sharp  leading  edge  is  specified  in  the  input  file,  rLE  is  set  equal  to 
0.01  in.  for  the  heating  computations. 

The  laminar  or  turbulent  status  of  the  flow  is  determined  by  computing  the 
Reynolds  number  based  on  the  leading  edge  diameter  and  the  free-stream  conditions. 


(4-67) 


If  ReD  is  below  some  user-specified  lower  limit,  ReL,  the  laminar  heat  transfer  is 
used.  If  ReD  is  greater  than  a  specified  upper  limit,  Reu,  the  turbulent  value  is  used. 
For  transitional  values  of  ReD  between  ReL  and  Reu,  the  heat  transfer  is  given  by 


(  \  . 

qw  “  qw,i  +  V  Rey  -  ReL  / 


(4-68) 


4.3  BODY 

For  points  on  the  body,  conditions  at  the  outer  edge  of  the  boundary  layer  must 
be  known.  In  the  case  of  pressure,  the  inviscid  surface  value  can  be  used  because,  to 
first  order,  the  pressure  gradient  across  the  boundary  layer  is  zero.  These  pressures 
are  known  from  the  appropriate  aeroprediction  methods.  During  the  course  of  the 
present  work,  some  refinements  were  made  to  the  blunt  body  pressure  methodology  of 
Reference  2.  These  changes  are  discussed  in  Appendix  A.  The  entropy  can  also  be 
determined  at  each  body  location.  If  the  bow  shock  is  attached  (sharp  body),  it  is 
simply  the  entropy  behind  the  bow  shock.  For  blunt  bodies,  an  iterative  boundary 
layer  mass  balance  technique  is  used  to  account  for  the  gradient  in  entropy  along  the 
boundary  layer  edge  created  by  the  curved  bow  shock.  This  approach  also  represents 
a  change  from  Reference  2  where  the  boundary  layer  edge  entropy  was  assumed  to  be 
constant  and  equal  to  the  entropy  behind  the  normal  shock  at  the  nose  tip.  This 
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constant  entropy  scenario  indeed  represents  the  true  inviscid  solution  because  the 
body  will  coincide  with  the  stagnation  streamline.  However,  the  streamlines  inter¬ 
secting  the  edge  of  a  boundary  layer,  which  is  increasing  in  thickness  in  the 
streamwise  direction,  originate  from  different  points  on  the  curved  bow  shock.  Thus, 
the  edge  entropy  varies  from  the  stagnation  value  at  the  nose  tip  to  some  limiting 
downstream  value  based  on  the  curvature  of  the  outer  bow  shock.  The  details  of  this 
entropy  computational  technique  are  described  in  Appendix  B. 

A  third  property,  the  total  enthalpy,  ho,  is  also  known  because  the  flow  is 
adiabatic.  The  remaining  fluid  properties  can  then  be  determined. 


4.3.1  Perfect  Gas 

As  before,  assume  a  constant  specific  heat,  then  proceed  as  follows: 


I 

c  = 

P  Y  -  1 

(4-69) 

exp(;r +  r +ln(ir)) 

(4-70) 

P  p  roo 

Pb 

Pb  ”b 

(4-71) 

hb='PTb 

(4-72) 

ab  =  ^  v  R  Tb 

(4-73) 

vb  =  V2(ho'V 

(4-74) 

M=  —  (4-75) 

b  ab 

zK  =  1  (4-76) 

b 
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4.3.2  Real  Gas 

Computations  proceed  along  the  same  lines  as  before  with  the  curve  fit  relations 
being  used. 


Pb  =  pb(pb-Sb) 

(4-77) 

hb  =  hb(pb>pb) 

(4-78) 

ab  =  ab(pb’sb) 

(4-79) 

Tb  =  Tb(pb*pb) 

(4-80) 

b  b  Pb 

(4-81) 

Vb  =  V2(ho“V 

(4-82) 

=  — 
b  ab 

(4-83) 

hb 

Yb=7 

eb 

(4-84) 

Pb 

zb_pbRTb 

(4-85) 

Prb  =  prb  (T1,-Pb) 

(4-86) 

pb  =  pb(Tb-pb) 

(4-87) 
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4.3.3  Eckert  Reference  Enthalpy  Flat  Plate  Heat  Transfer 

Given  the  appropriate  fluid  properties,  the  laminar  and  turbulent  heat  transfer 
rates  are  computed  asl3 


q  ,  =  0.332  (Pr*) 

W,1 


P  vi. 
-0.667  b 


N. 


and 


qwt  =  0.185  (Pr*) 


-0.667 


P 


m)2 


Re_\\2.584 

The  *  conditions  are  evaluated  at  the  reference  enthalpy,  h*. 


h  =  0.5  (h  +  h  )  +  0.22  (h  -h  ) 

we  aw  e 


For  a  perfect  gas,  the  reference  temperature,  T*,  may  be  used. 


T  =  0.5  (T  +  T  )  +  0.22  (T  -T  ) 

we  aw  e 


(4-88) 


(4-89) 


(4-90) 


(4-91) 


The  e  subscript  once  again  denotes  the  conditions  at  the  edge  of  the  boundary 
layer  and,  in  this  case,  corresponds  to  the  b  subscript.  Ni  and  Nt  are  the  laminar  and 
turbulent  Mangier  transformation  factors  that  are  equal  to  3.0  and  2.0,  respectively. 
These  factors  provide  an  approximation  to  three-dimensional  results  even  though  the 
basic  heat  transfer  relations  were  originally  derived  for  a  flat  plate.  The  degree  of 
turbulence  is  determined  by  first  computing  a  local  Reynolds  number  based  on  the 
boundary  layer  running  length,  s,  and  the  local  reference  conditions. 


Re* 


p*yb8 

* 

P 


(4-92) 


As  before,  if  Re*  is  less  than  a  user-specified  lower  limit,  ReL,  the  flow  is 
assumed  to  be  fully  laminar.  For  Re*  values  greater  than  a  specified  upper  limit, 
Reu,  the  flow  is  assumed  to  be  fully  turbulent.  For  intermediate  transitional  values 
of  Re*  between  ReL  and  Reu,  the  heat  transfer  is  given  by 


__  (Be  -  ReL  ^ 
qw  qw.l  +  l  Re..  -  Re. 


Re  —  Re 

L 


(4-93) 
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The  manner  in  which  s  is  determined  within  the  program  is  strictly  correct  only 
for  zero  angle  of  attack.  It  is  computed  along  constant  circumferential  lines  and  is 
equivalent  to  the  streamwise  body  coordinate,  referenced  to  the  stagnation  point  as 
origin.  In  actuality,  s  should  be  measured  along  surface  streamlines.  At  zero  angle  of 
attack,  the  streamlines  coincide  with  the  body  coordinate  system,  but  as  angle  of 
attack  increases,  the  windward-side  streamlines  tend  to  sweep  around  the  body 
toward  the  leeward  side.  To  include  this  effect  in  the  aeroprediction  heating  analysis 
would  have  involved  a  considerable  penalty  in  computational  time,  so  the  approxi¬ 
mate  approach  was  chosen.  At  high  angles  of  attack,  this  approximation  can 
adversely  impact  the  results. 


4.4  CONTROL  SURFACES 

The  control  surface  heating  is  determined  in  a  manner  exactly  analagous  to 
that  for  the  body.  The  boundary  layer  edge  conditions  are  available  from  the 
appropriate  aeroprediction  subroutines  and  are  incorporated  directly  into  the  Eckert 
Reference  Enthalpy  Flat  Plate  relations  given  by  Equations  4-88  and  4-89.  The 
Mangier  transformation  factors,  Ni  and  Nt,  are  both  set  to  one,  since  the  flow  is 
essentially  two-dimensional  along  the  surface  of  fins  or  wings  and  the  original  flat 
plate  form  of  the  equations  is  appropriate.  The  reference  properties  are  computed  as 
for  the  body  surface,  and  the  laminar/turbulent  character  of  the  flow  is  determined  by 
comparing  the  local  reference  Reynolds  number  to  the  specified  upper  and  lower 
limits.  If  the  flow  is  transitional,  Equation  4-93  is  used.  The  boundary  coordinate,  s, 
is  measured  straight  back  at  a  constant  radial  distance  from  the  leading  edge. 
Because  there  can  be  some  inboard  or  outboard  component  to  the  flow,  this  is  only  an 
approximation  to  the  actual  streamline  distance.  The  error  involved  in  this  simplifi¬ 
cation  is  usually  quite  small,  and  any  attempt  to  compensate  for  it  would  not  be  cost 
effective  for  present  purposes. 
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5.0  RESULTS 


Validation  of  the  new  heat  transfer  methodology  was  carried  out  in  two  steps. 
The  first  involved  comparing  results  from  aeroprediction  with  those  from  the 
MINIVER  code  for  stagnation  point  and  stagnation  line  flows  and  also  for  sharp  cone 
and  sharp  wedge  flows  at  zero  angle  of  attack.  MINIVER  models  these  specialized 
configurations  quite  well,  providing  a  good  check  for  both  heat  transfer  and  inviscid 
boundary  layer  edge  conditions.  The  latter  comparisons  served  as  further  validation 
of  the  real-gas  model  previously  incorporated  into  the  aeroprediction  code.  In  the 
second  step  of  the  validation  process,  computations  were  performed  for  spherically 
blunted  cones  at  angle  of  attack.  The  resulting  heat  transfer  rates  were  compared  to 
values  generated  by  more  sophisticated  numerical  methods  and  to  experimental 
measurements. 


5.1  STAGNATION  POINT 

The  free-stream  conditions  for  these  computations  correspond  to  an  altitude  of 
82,500  ft  and  Mach  numbers  of  7.73, 10,  and  15.  The  nose  radius  was  0.01  in.  and  the 
surface  temperature  was  set  to  580  °R.  Aeroprediction  and  MINIVER  results  for 
boundary  layer  edge  pressure,  density,  temperature,  and  enthalpy,  as  well  as 
recovery  temperature  and  heat  transfer  rate,  are  shown  in  Tables  5-1  through  5-3.  In 
general,  results  from  the  two  codes  differ  by  less  than  2  percent  for  all  variables.  The 
only  exceptions  are  the  heat  transfer  rates  at  Mach  10  and  15  with  differences  of 
2.4  and  3.6  percent,  respectively.  The  primary  reason  for  these  discrepancies  lies  in 
the  differing  values  of  viscosity  coefficient.  MINIVER  uses  a  modified  Sutherland 
Law  formulation  to  compute  the  coefficient  of  viscosity,  while  aeroprediction  uses 
curve  fits.8,9  As  the  temperature  increases,  these  two  approaches  give  increasingly 
divergent  values  with  the  difference  approaching  33  percent  at  Mach  15.  The 
magnitude  of  this  variation  is  not  fully  reflected  in  the  heating  rates  because  of  their 
square  root  dependency  on  viscosity  coefficient.  There  are  also  compensating 
differences  in  the  Prandtl  number.  MINIVER  uses  a  constant  value  of  0.72  while 
aeroprediction  once  again  relies  on  the  curve  fits.  Because  the  curve  fits  have  a 
stated  range  of  validity  of  up  to  25,000  °R,  they  should  offer  the  more  accurate 
approach. 
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TABLE  5-1.  NOSE  STAGNATION  POINT  COMPARISONS,  M*,  =  7.73 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

4080 

4120 

Density  (slug/ft3) 

0.0005335 

0.000540 

Temperature  (°R) 

4501 

4420 

Enthalpy  (ft2/sec2) 

31121047 

31038431 

Recovery  temperature  (°R) 

4501 

4419 

Heat  transfer  (BTU/ft2-sec) 

2170 

2202 

TABLE  5-2.  NOSE  STAGNATION  COMPARISONS,  M»  =  10 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

6894 

6900 

Density  (slug/ft3) 

0.0006453 

0.000658 

Temperature  (°R) 

5877 

5924 

Enthalpy  (ft2/sec2) 

50480874 

50311629 

Recovery  temperature  (°R) 

5877 

5924 

Heat  transfer  (BTU/ft2-sec) 

4856 

4974 

TABLE  5-3.  NOSE  STAGNATION  COMPARISONS,  M.  =  15 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

15640 

15800 

Density  (slug/ft3) 

0.0008021 

0.000816 

Temperature  (°R) 

9501 

9453 

Enthalpy  (ft2/sec2) 

110577829 

110214264 

Recovery  temperature  (°R) 

9501 

9453 

Heat  transfer  (BTU/ft2-sec) 

17254 

17900 

5-2 


NSWCDD/TR-93/29 


5.2  STAGNATION  LINE 

The  second  configuration  considered  was  the  stagnation  line  of  a  swept  cylinder. 
As  in  the  first  case,  free-stream  conditions  correspond  to  an  altitude  of  82,500  ft  at 
Mach  numbers  of  7.73,  10,  and  15.  The  cylinder  radius  was  0.01  in.  with  a  surface 
temperature  of  580  °R.  The  sweep  angle  was  53.13  deg.  Representative  results  from 
these  runs  are  shown  in  Tables  5-4  through  5-6.  Agreement  is  good  overall,  with 
maximum  differences  of  2.5  percent  or  less  except  for  heat  transfer  rates  where 
discrepancies  run  somewhat  higher  for  the  reasons  discussed  previously.  For  the  two 
higher  Mach  numbers,  the  laminar  heating  rates  shown  are  higher  than  the 
corresponding  turbulent  rates.  This  situation  can  arise  when  the  turbulent  relations 
are  applied  at  a  Reynolds  number  that  is  below  their  range  of  validity.  The  results 
were  still  included  in  the  tables  because  they  provide  further  confirmation  that  the 
heating  rate  computations  in  aeroprediction  are  behaving  as  they  should. 


TABLE  5-4.  SWEPT  CYLINDER  STAGNATION  LINE  COMPARISONS,  =  7.73 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

1470 

1480 

Density  (slug/ft3) 

0.0004266 

0.00430 

Temperature  (°R) 

1992 

2004 

Enthalpy  (ft2/sec2) 

12745350 

12704919 

Recovery  temperature  (°R) 

4013 

4025 

Laminar  heat  transfer  (BTU/ft2-sec) 

492 

480 

Turbulent  heat  transfer  (BTU/ft2-sec) 

830 

802 

TABLE  5-5.  SWEPT  CYLINDER  STAGNATION  LINE  COMPARISONS,  M„  =  10 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

2453 

2490 

Density  (slug/ft3) 

0.0004742 

0.000486 

Temperature  (°R) 

2978 

2978 

Enthalpy  (ft2/sec2) 

19711317 

19652980 

Recovery  temperature  (°R) 

5433 

5516 

Laminar  heat  transfer  (BTU/ft2-sec) 

1870 

1817 

Turbulent  heat  transfer  (BTU/ft2-sec) 

1125 

1200 
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TABLE  5-6.  SWEPT  CYLINDER  STAGNATION  LINE  COMPARISONS,  M»  =  15 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

5560 

5570 

Density  (slug/ft3) 

0.0005942 

0.000606 

Temperature  (°R) 

5286 

5286 

Enthalpy  (ft2/sec2) 

41346321 

41235184 

Recovery  temperature  (°R) 

8701 

8482 

Laminar  heat  transfer  (BTU/ft2-sec) 

6727 

6550 

Turbulent  heat  transfer  (BTU/ft2-sec) 

4526 

4700 

As  a  further  check,  two  additional  cylinder  radii  of  0.5  and  1.0  in.  were  run. 
Heat  transfer  rates  for  these  cases  are  shown  in  Tables  5-7  and  5-8.  Agreement 
between  the  two  methods  is  consistent  with  the  previous  results. 

TABLE  5-7.  SWEPT  CYLINDER  STAGNATION  LINE 
COMPARISONS,  RADIUS  =  0.5  IN. 


Man 

Code 

Laminar 

Heat  Transfer 
(BTU/ft2-sec) 

Turbulent 

Heat  Transfer 
(BTU/ft2-sec) 

7.73 

Aeroprediction 

117 

MINIVER 

113 

hBHI 

10.0 

Aeroprediction 

264 

514 

MINIVER 

257 

550 

15.0 

Aeroprediction 

951 

■■■ 

MINIVER 

926 

TABLE  5-8.  SWEPT  CYLINDER  STAGNATION  LINE 
COMPARISONS,  RADIUS  =  1.0  IN. 


Mo, 

Code 

Laminar 

Heat  Transfer 
(BTU/ft2-sec) 

Turbulent 

Heat  Transfer 
(BTU/ft2-sec) 

7.73 

Aeroprediction 

83 

180 

MINIVER 

80 

190 

10.0 

Aeroprediction 

187 

MINIVER 

182 

15.0 

Aeroprediction 

673 

1802 

MINIVER 

650 

1900 
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5.3  SHARP  CONE 

For  the  sharp-cone  computations,  free-stream  conditions  were  again  taken  to  be 
Mach  numbers  of  7.73, 10,  and  15  at  an  altitude  of  82,500  ft.  The  cone  half  angle  was 
30  deg,  and  the  heat  transfer  was  computed  at  a  point  1.75  in.  from  the  nose  tip 
stagnation  point.  Surface  temperatures  were  836  °R  at  Mach  7.73,  1200  °R  at 
Mach  10,  and  3658  °R  at  Mach  15.  Representative  results  from  aeroprediction  and 
MINIVER  for  these  conditions  are  shown  in  Tables  5-9  through  5-11.  Agreement  is 
good,  with  differences  being  well  under  2  percent  in  most  instances.  The  exception  is 
the  heat  transfer  rate  at  the  highest  Mach  number  where  the  turbulent  values  differ 
by  4.6  percent.  The  other  heat  transfer  rates  closely  agree,  a  situation  brought  about 
by  the  lower  temperatures  involved  in  these  computations  relative  to  the  stagnation 
conditions  and  a  corresponding  closer  agreement  in  viscosity  coefficients. 


TABLE  5-9.  30-DEG  SHARP-CONE  COMPARISONS,  Moo  =  7.73 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

1201 

1210 

Density  (slug/ft3) 

0.0003888 

0.000396 

Temperature  (°R) 

1787 

1782 

Enthalpy  (ft2/sec2) 

11332032 

11169005 

Recovery  temperature  (°R) 

3977 

3986 

Laminar  heat  transfer  (BTU/ft2-sec) 

80.4 

81.0 

Turbulent  heat  transfer  (BTU/ft2-sec) 

224 

221 

TABLE  5-10.  30-DEG  SHARP-CONE  COMPARISONS,  Moo  =  10 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

1957 

1980 

Density  (slug/ft3) 

0.0004384 

0.000450 

Temperature  (°R) 

2574 

2561 

Enthalpy  (ft2/sec2) 

16806866 

16636016 

Recovery  temperature  (°R) 

5367 

5460 

Laminar  heat  transfer  (BTU/ft2-sec) 

181 

180 

Turbulent  heat  transfer  (BTU/ft2-sec) 

511 

510 
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TABLE  5-11.  30-DEG  SHARP-CONE  COMPARISONS,  M«  =  15 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

4292 

4350 

Density  (slug/ft3) 

0.0005449 

0.000547 

Temperature  (°R) 

4608 

4606 

Enthalpy  (ft2/sec2) 

32248110 

33104186 

Recovery  temperature  (°R) 

8576 

8387 

Laminar  heat  transfer  (BTU/ft2-sec) 

571 

565 

Turbulent  heat  transfer  (BTU/ft2-sec) 

1652 

1579 

5.4  SHARP  WEDGES 

Once  again,  computations  were  performed  for  Mach  numbers  of  7.73, 10,  and  15 
at  an  altitude  of  82,500  ft.  The  wedge  half  angle  was  4  deg,  and  the  wall  temperature 
was  a  constant  580  °R.  Heating  rates  were  computed  at  a  point  0.655  in.  from  the 
stagnation  line.  Results  are  shown  in  Tables  5-12  through  5-14.  Agreement  is  once 
again  good  with  boundary  layer  edge  conditions  differing  by  less  than  3  percent  and 
heat  transfer  rates  by  4  percent  or  less. 


TABLE  5-12.  4-DEG  SHARP-WEDGE  COMPARISONS,  Mx  =  7.73 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

107 

109 

Density  (slug/ft3) 

0.0001256 

0.000127 

Temperature  (°R) 

496 

500 

Enthalpy  (ft2/sec2) 

2982836 

2992163 

Recovery  temperature  (°R) 

3763 

3810 

Laminar  heat  transfer  (BTU/ft2-sec) 

26.6 

25.6 

Turbulent  heat  transfer  (BTU/ft2-sec) 

53.4 

51.0 

NSWCDD/TR-93/29 


TABLE  5-13.  4-DEG  SHARP-WEDGE  COMPARISONS,  M®  =  10 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

130 

134 

Density  (slug/ft3) 

0.0001424 

0.000146 

Temperature  (°R) 

530 

536 

Enthalpy  (ft2/sec2) 

3184599 

3211866 

Recovery  temperature  (°R) 

4858 

4967 

Laminar  heat  transfer  (BTU/ft2-sec) 

53.2 

51.5 

Turbulent  heat  transfer  (BTU/ft2-sec) 

104 

100 

TABLE  5-14.  4-DEG  SHARP-WEDGE  COMPARISONS,  M*,  =  15 


Parameter 

Aeroprediction 

MINIVER 

Pressure  (lb/ft2) 

203 

208 

Density  (slug/ft3) 

0.0001849 

0.000189 

Temperature  (°R) 

639 

641 

Enthalpy  (ft2/sec2) 

3707810 

3842665 

Recovery  temperature  (°R) 

7303 

7340 

Laminar  heat  transfer  (BTU/ft2-sec) 

162 

158 

Turbulent  heat  transfer  (BTU/ft2-sec) 

301 

300 

5.5  SPHERICALLY  BLUNTED  CONES 

The  first  configuration  considered  consisted  of  a  15-deg  half-angle  cone  with  a 
nose  radius  of  0.375  in.  Free-stream  pressure  and  temperature  were  2.66  lb/ft2  and 
89.971  °R.  The  Mach  number  was  10.6  at  an  angle  of  attack  of  20  deg.  The  wall 
temperature  was  540  °R.  Heat  transfer  rates  along  the  windward  symmetry  plane 
are  shown  as  a  function  of  distance  along  the  body  in  Figure  5-1.  For  comparison 
purposes,  numerical  results  are  shown  from  Reference  14  and  corresponding 
experimental  results  from  Reference  15.  The  methodology  used  in  Reference  14  was 
much  more  involved  than  the  approach  employed  in  aeroprediction  with  the  inclusion 
of  complex  streamline  tracking  methods  and  the  axisymmetric  analogue  technique  to 
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model  streamline  spreading  effects.  Aeroprediction  and  MINIVER  tend  to  under¬ 
predict  the  experimental  data  and  the  computational  results  of  Reference  10  by  15  to 
20  percent,  a  performance  that  is  quite  creditable  considering  the  simplified  solution 
approach.  At  the  nose  tip  stagnation  point,  aeroprediction  gives  a  heat  transfer  rate 
of  33.99  BTU/ft2-sec  while  Reference  14  predicts  35.94  BTU/ft2-sec,  which  again  is 
good  agreement. 


HEAT 

TRANSFER 

RATE 

(BTU/ft2-sec) 


x 


FIGURE  5-1.  HEAT  TRANSFER  RATES  FOR  0.375-IN.  NOSE  RADIUS, 
15-DEG  HALF-ANGLE  CONE  AT  a  =  20  DEG 
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The  second  configuration  considered  was  a  spherically  blunted,  15-deg  half¬ 
angle  cone  with  a  nose  radius  of  1.1  in.  Free-stream  pressure  and  temperature  were 
2.66  lb/ft2  and  89.971  °R,  respectively.  The  Mach  number  was  10.6  and  angles  of 
attack  of  5  and  10  deg  were  run.  The  wall  temperature  was  a  constant  540  °R. 
Figures  5-2  and  5-3  show  heating  rate  predictions  on  the  windward  symmetry  plane 
as  a  function  of  position  along  the  body  for  the  5-  and  10-deg  angle-of-attack  cases, 
respectively.  Also  shown  are  results  from  AEROHEAT,  an  engineering  type  code 
that  uses  streamline  tracking  and  axisymmetric  analogue  methodology,  and  from  a 
second  approximate  code  that  uses  more  sophisticated  streamline  tracking,  as  well  as 
a  more  exact  boundary  layer  formulation.  16  MINIVER  results  are  also  plotted  along 
with  experimental  data.15  Aeroprediction  does  well  beyond  about  5  nose  radii 
downstream  but  tends  to  overpredict  heating  rates  near  the  nose  tip.  The  region  of 
overprediction  is  somewhat  less  at  the  higher  angle  of  attack.  This  discrepancy  is 
believed  to  be  caused  by  the  use  of  heating  rate  expressions  that  are  based  on  flat 
plate  relations  and  thus  are  incapable  of  modeling  the  effects  of  the  rapidly  changing 
boundary  layer  edge  conditions  in  this  region.  The  aeroprediction  results  are  slightly 
better  than  those  from  MINIVER  because  of  the  inclusion  of  variable  entropy  effects. 
Figures  5-4  through  5-7  show  the  computed  variation  of  heating  rates  in  the 
circumferential  direction  around  the  body.  The  <J>  =  0  location  corresponds  to  the 
windward  symmetry  plane.  Figures  5-4  and  5-5  illustrate  results  for  a  5-deg  angle  of 
attack  at  points  4.86  and  10.13  nose  radii  downstream  of  the  nose  tip,  respectively. 
At  4.86  nose  radii,  the  aeroprediction  results  are  uniformly  high,  a  result  that  could 
be  anticipated  because  of  the  tendency  to  overpredict  heating  rates  in  the  nose  region. 
At  10.13  nose  radii,  agreement  with  both  numerical  and  experimental  results  is  good, 
with  the  greater  differences  occurring  on  the  leeward  side  of  the  body.  Figures  5-6 
and  5-7  show  the  10-deg  angle-of-attack  results  for  the  same  streamwise  locations. 
At  the  4.86  nose  radii  station,  except  for  the  leeward  plane,  the  aeroprediction  results 
are  actually  more  in  agreement  with  the  experiment  than  are  the  predictions  of  the 
more  sophisticated  numerical  methodology  of  Reference  16.  At  10.13  nose  radii,  the 
agreement  is  not  as  good  but  is  still  within  15  percent  of  both  experimental  data  and 
other  computational  predictions.  No  results  are  shown  from  MINIVER  for  these 
circumferentially  varying  cases  because  it  does  not  explicitly  provide  a  capability  to 
handle  three-dimensional  problems. 
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FIGURE  5-2.  HEAT  TRANSFER  RATES  FOR  1.1-IN.  NOSE  RADIUS, 
15-DEG  HALF-ANGLE  CONE  AT  a  =  5  DEG 
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FIGURE  5-3.  HEAT  TRANSFER  RATES  FOR  1.1-IN.  NOSE  RADIUS, 
15-DEG  HALF-ANGLE  CONE  AT  a  =  10  DEG 
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FIGURE  5-4.  CIRCUMFERENTIAL  VARIATION  OF  HEATING  RATES  ON  1.1-IN. 
NOSE  RADIUS,  15-DEG  HALF-ANGLE  CONE  AT  a  =  5  DEG,  xJr^  =  4.86 
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FIGURE  5-5.  CIRCUMFERENTIAL  VARIATION  OF  HEATING  RATES  ON  1.1-IN. 
NOSE  RADIUS,  15-DEG  HALF-ANGLE  CONE  AT  a  =  5  DEG,  x/r*  =  10.13 
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FIGURE  5-6.  CIRCUMFERENTIAL  VARIATION  OF  HEATING  RATES  ON  1.1-IN. 
NOSE  RADIUS,  15-DEG  HALF-ANGLE  CONE  AT  a  =  10  DEG,  x/i*  =  4.86 
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FIGURE  5-7.  CIRCUMFERENTIAL  VARIATION  OF  HEATING  RATES  ON  1.1-IN. 
NOSE  RADIUS,  15-DEG  HALF-ANGLE  CONE  AT  a  =  10  DEG,  x/rn  =  10.13 
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6.0  SUMMARY 


Methodology  has  been  incorporated  into  the  aeroprediction  code  to  permit 
computation  of  heat  transfer  rates  and  recovery  temperatures  as  part  of  the  high 
Mach  number  solution.  The  aeroheating  techniques  implemented  include  a 
simplified  Fay-Riddell  model  for  the  nose  stagnation  point,  the  Beckwith-Gallagher 
swept-cylinder  method  (modified  to  include  real-gas  effects)  for  control  surface 
leading  edges,  and  the  Eckert  reference  enthalpy  flat  plate  method  for  all  other 
points  on  the  body  and  control  surfaces.  A  mass  balance  technique  has  been  added  to 
determine  the  appropriate  entropy  value  to  use  as  a  boundary  layer  edge  condition 
when  solutions  are  computed  for  blunt  bodies. 

During  the  course  of  this  work,  refinements  were  made  to  the  aeroprediction 
pressure  methodology  to  achieve  a  better  match  point  between  the  Improved 
Modified  Newtonian  Theory  used  on  the  forward  portion  of  blunt  bodies  and  the 
Second-Order  Shock-Expansion  Theory  solution  downstream  of  the  nose  region. 
Further  modifications  were  made  to  remove  discontinuities  in  the  circumferential 
pressure  distribution  resulting  in  improvements  in  overall  normal  and  axial  force 
coefficients,  especially  for  blunted  cones. 

Comparisons  of  heat  transfer  predictions  from  the  updated  code  with  those  from 
more  sophisticated  engineering  analysis  methods  and  with  experimental  data  show, 
in  general,  good  agreement.  There  is  a  tendency  to  overpredict  the  heat  transfer  near 
the  nose  tip— -a  situation  most  likely  arising  from  the  inability  of  the  methods  used  to 
account  for  the  rapidly  changing  boundary  layer  edge  conditions  in  this  region. 

Previously,  most  aeroheating  analysis  at  NSWCDD  has  been  done  using  either 
the  MINIVER  code  or  a  combination  ZEUS/boundary  layer  (ZEUSBL)  approach. 
MINIVER  has  the  advantage  of  being  easy  to  implement  and  fast  to  run,  producing 
time  histories  of  surface  temperatures  at  selected  points  with  minimal  effort.  It 
suffers  from  an  inability  to  model  three-dimensional  effects  in  an  explicit  manner  or 
to  account  for  a  surface  entropy  layer.  ZEUSBL  can  model  three-dimensional  effects 
well,  but  at  the  expense  of  much  greater  computational  effort.  It  does,  however,  have 
the  drawback  of  lacking  a  variable  entropy  model  in  its  present  form.  The  present 
method  offers  an  approach  to  the  aeroheating  problem  that,  while  not  as 
sophisticated  or  rigorous  as  ZEUSBL,  represents  a  definite  improvement  over 
MINIVER  in  analyzing  bodies  at  angle  of  attack.  It  is  also  the  superior  method  if 
surface  entropy  layer  effects  are  significant.  This  added  capability  is  gained  at  the 
expense  of  only  a  modest  increase  in  computational  effort. 
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8.0  SYMBOLS 

Speed  of  sound  (ft/sec) 

Axial  force  coefficient 

Normal  force  coefficient 

Constant  pressure  specific  heat  (ft2/sec2-°R) 

Diameter  (ft) 

Internal  energy  (ft2/sec2) 

Temperature-based  heat  transfer  coefficient  (ft-lb/ft2-sec-°R) 
Enthalpy  (ft2/sec2) 

Mach  number 

Mangier  transformation  factor 
Prandtl  number 
Pressure  (lb/ft2) 

Heat  transfer  (ft-lb/ft2-sec  or  BTU/ft2-sec) 

Gas  constant  (ft2/sec2-°R) 

Shock  radius  of  curvature  at  vertex  (ft) 

Reynolds  number 
Body  or  shock  radius  (ft) 

Recovery  factor 
Entropy  (ft-lb/slug-°R) 
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SUBSCRIPTS 

aw 
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Boundary  layer  running  length  (ft) 

Temperature  (°R) 

Velocity  in  streamwise  direction  (ft/sec) 

Velocity  (fit/sec) 

Streamwise  coordinate  (ft) 

Compressibility  factor 
Angle  of  attack  (deg) 

Limiting  shock  wave  angle  (deg) 

Specific  heat  ratio 
Boundary  layer  thickness  (ft) 

Boundary  layer  displacement  thickness  (ft) 

Cone  half  angle  (deg) 

Enthalpy -based  heat  transfer  coefficient  (slug/ft2-sec) 
Local  body  slope  in  streamwise  direction  (deg) 

Control  surface  leading  edge  sweep  (deg) 

Viscosity  coefficient  (slug/ft-sec) 

Density  (slug/ft3) 

Circumferential  angular  coordinate  (deg) 

Shock  standoff  distance  (ft) 


Adiabatic  wall 

Body 

Diameter 
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Downstream 
Boundary  layer  edge 
Lower  transition  limit 
Leading  edge 
Laminar 
Match  point 
Nose  tip 

Normal  direction 
Parallel 

Shock  conditions 

Based  on  boundary  layer  running  length 
Turbulent 

Upper  transition  limit 
Wall  or  surface 
Stagnation  conditions 
Downstream  of  shock  or  expansion 
Free-stream  conditions 


SUPERSCRIPTS 

'  Inviscid  surface  conditions  behind  parallel  shock 

*  Eckert  reference  conditions 


8-3 


NSWCDD/TR-93/29 


APPENDIX  A 

MODIFICATIONS  TO  PRESSURE  METHODOLOGY 
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For  blunt  bodies,  the  pressure  prediction  methodology  of  Reference  A-l  involves 
a  two-step  procedure  in  which  an  improved  modified  Newtonian  Technique  (IMNT)  is 
applied  to  the  nose  region  and  Second-Order  Shock-Expansion  Theory  (SOSET)  is 
used  for  the  remainder  of  the  body.  Originally,  the  match  point  between  the  IMNT 
and  SOSET  regions  was  chosen  to  be  the  location  at  which  the  local  body  slope  had 
decreased  to  25.95  deg.  It  was  found  that  using  this  fixed  match  point  often  resulted 
in  pressure  discontinuities  of  10  to  15  percent.  To  eliminate  this  inconsistency,  a 
parametric  study  was  undertaken  to  determine  the  actual  slope  at  which  IMNT  and 
SOSET  produced  equivalent  pressures  over  a  range  of  Mach  numbers.  The  resulting 
Mach  number/match  point  correlations  were  fit  to  a  quadratic  expression  for  match 
point  as  a  function  of  Mach  number. 

0M  =  -0.028422619  M2  +  1.028273808  M  +  13.35357144 

Determining  the  match  point  in  this  manner  decreases  the  pressure  dis¬ 
continuity  to  less  than  2  percent  in  the  Mach  number  range  from  6  to  20,  as  shown  in 
Table  A-l.  The  error  will  be  greater  below  Mach  6,  but  this  pressure  methodology 
would  not  normally  be  used  for  the  lower  Mach  numbers.  The  results  shown  in 
Table  A-l  represent  zero  angle-of-attack  values.  In  aeroprediction,  the  zero  incidence 
solution  is  computed  first  and  then  angle-of-attack  effects  are  added  by  adjusting 
these  pressures  with  an  approximate  analytical  relation  involving  angle-of-attack, 
local  body  slope  and  body  circumferential  location.  During  this  transformation 
process,  the  match  point  errors  can  sometimes  increase  slightly,  but  rarely  do  they 
exceed  2  percent,  and  in  no  case  were  they  found  to  exceed  3  percent. 

A  second  problem  with  pressure  discontinuities  was  encountered  that  was 
associated  with  the  manner  in  which  the  angle-of-attack  corrections  mentioned  were 
applied.  In  Reference  A-l,  it  was  found  that  the  use  of  a  modified  form  of  the  pressure 
adjustment  relations  gave  better  results  in  the  leeward  planes  for  nonconical  bodies. 
Unfortunately,  it  was  found  that  using  this  modified  formulation  on  the  leeward  side 
of  blunted  cones  resulted  in  sharp  pressure  discontinuities  in  the  circumferential 
direction.  The  code  was  modified  to  apply  the  pressure  correction  uniformly  at  all 
points  for  conical  bodies  and  to  slightly  extend  the  region  of  modified  corrections  for 
the  nonconical  bodies  from  30  to  60  deg  away  from  the  leeward  body  centerline.  This 
procedure  eliminated  the  discontinuities  and,  in  terms  of  total  coefficients  as  shown 
in  Table  A-2,  gave  much  better  results  for  cones  and  slightly  better  ones  for  ogives. 


A-l  Moore,  F.  G.;  Armistead,  M.  A.;  Rowles,  S.  H.;  and  DeJamette,  F.  R..  Second-Order  Shock- 
Expansion  Theory  Extended  to  Include  Real  Gas  Effects ,  NAVSWC  TR  90-683,  Feb  1992, 
NAVSWC,  Dahlgren,  VA. 
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TABLE  A-l.  MATCH  POINT  VARIATION 
WITH  MACH  NUMBER 


Mach  Number 

Body  Slope  at 
Match  Point 
(deg) 

Pressure  Error 
(percent) 

6 

18.5 

0.44 

8 

20.1 

1.75 

10 

21.0 

1.38 

12 

21.6 

0.52 

15 

22.0 

0.62 

18 

22.4 

0.31 

20 

22.6 

0.51 

TABLE  A-2.  EFFECT  OF  PRESSURE  MODIFICATIONS  AT 
ANGLE  OF  ATTACK-M  =  15,  ANGLE  OF  ATTACK  =  0  DEG 

20-Percent  Blunt  20-Percent  Blunt 

15-Deg  Cone  Von  Karman  Ogive 


Cn  I  CA  I  CN 


Reference  A-l 

0.257 

0.219 

0.314 

0.149 

Present  method 

0.313 

0.201 

0.308 

0.150 

ZEUS 

0.297 

0.190 

0.284 

0.112 
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To  apply  any  boundary  layer  solution  technique,  it  is  first  necessary  to 
determine  appropriate  outer  edge  boundary  conditions.  One  method  of  accomplishing 
this  task  is  to  first  obtain  an  inviscid  solution  for  the  flow  about  the  body  of  interest 
and  then  apply  the  surface  conditions  from  this  model  as  boundary  layer  edge 
conditions.  For  pointed  bodies  with  attached  shocks  of  uniform  slope,  this  method 
provides  physically  consistent  solutions.  However,  in  the  case  of  blunted  objects,  the 
situauon  becomes  more  involved.  The  true  inviscid  surface  solution  implies  that  the 
entropy-  along  the  body  surface  is  a  constant  equal  to  the  value  behind  the  normal 
shock  at  the  forward  stagnation  point.  This  situation  results  because  of  the  con¬ 
stancy  of  entropy  along  streamlines  and  the  fact  that  the  body  surface  is  coincident 
with  the  stagnation  streamline.  Two  factors  make  the  use  of  the  constant  normal 
shock  entropy  as  the  boundary  layer  edge  condition  inappropriate.  First,  the  bow 
shock  ahead  of  a  blunt  object  is  highly  curved  in  the  vicinity  of  the  body.  Thus,  the 
streamlines  lying  near  the  body  surface  cross  the  shock  at  locations  with  different 
slopes  resulting  in  different  entropy  values  along  each  streamline.  An  entropy  layer 
is  created  that  flows  back  over  the  surface.  The  second  complicating  factor  is  the 
finite  thickness  of  the  boundary  layer  that  increases  in  the  streamwise  direction.  It  is 
the  value  of  entropy  along  the  streamline  that  intersects  the  outer  boundary  layer 
that  determines  the  appropriate  edge  properties  to  use  in  the  viscous  solution.  Only 
if  the  boundary  layer  were  of  zero  thickness  would  the  edge  entropy  be  the  normal 
shock  value.  This  situation  is  illustrated  in  Figure  B-l. 


FIGURE  B-l.  DETERMINATION  OF  BOUNDARY  LAYER  EDGE  ENTROPY 
ON  BLUNT  BODIES  WITH  HIGHLY  CURVED  BOW  SHOCKS 
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To  determine  this  local  entropy  value,  a  mass  balance  technique  can  be  applied. 
The  mass  flux  through  a  basically  trapezoidal  element  bounded  by  the  body  surface 
and  the  boundary  layer  edge  with  a  circumferential  angular  displacement  is 

1  (B-l) 

m,  =  -  p  V  (r.  +  8*  +  8)  (8  -  8*)  A<J> 
b  2  c  e  ° 


where 


pe  =  Density  at  boundary  layer  edge 

Ve  =  Velocity  at  boundary  layer  edge 

rb  =  Local  body  radius 

8  =  Local  boundary  layer  thickness 

8*  =  Local  boundary  layer  displacement  thickness 

At  the  shock,  the  mass  flux  through  a  triangular  element  A<J>  in  circumferential 
displacement  and  rs  in  radius  is 


(B-2) 


where 


poo  —  Free-stream  density 
V  00  ■  Free-stream  density 

If  mb  and  ms  are  equated,  it  is  possible  to  solve  for  rs,  which  will  be  the  radial 
location  on  the  shock  of  the  streamline  that  intersects  the  boundary  layer  at  the  local 
body  point.  This  gives 


r 


s 


(2  rb  +  8*  +  8)  (8  -  8*) 


(B-3) 


The  determination  of  rs  is  an  iterative  process,  because  the  local  boundary  lryer 
thickness  is  dependent  on  the  local  entropy.  A  value  for  entropy  must  be  assumed, 
and  in  conjunction  with  the  local  pressure  value  that  is  known  and  is  not  affected  by 
the  boundary  layer  (to  a  first  approximation),  the  other  local  thermodynamic 
properties  can  be  determined.  A  boundary  layer  thickness  is  then  computed.  At  the 
present  time,  this  computation  is  done  by  basic  incompressible  relations  with  the  use 
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of  Eckert’s  reference  enthalpy  properties  to  account  for  variations  across  the 
boundary  layer.  For  laminar  flow 


5.0  s 

*l“  VrT 

a 


(B-4) 


8 


1 


1.729  s 


(B-5) 


where  s  is  the  local  boundary  layer  coordinate  and  Res  is  the  Reynolds  number  based 
on  s  and  the  local  reference  enthalpy  properties.  For  turbulent  flow 


0.37  b 
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For  local  Reynolds  numbers  below  some  lower  value  ReL,  the  laminar  values  are 
used.  For  local  Reynolds  numbers  above  an  upper  limit  Reu,  the  turbulent  value  is 
used.  For  transitional  cases  with  ReL<Res<Reu>  the  boundary  layer  thickness  is 
given  by 


H  _R®L 

8  =  8.+  - - —  (8.  -  8.) 

1  Re  —  Re  t  1 


(B-8) 


The  Reynolds  number  adjusted  value  of  8*  is  determined  in  a  similar  fashion. 
Knowing  8  and  8*  for  the  assumed  value  of  entropy,  rs  can  be  determined.  The  shock 
slope  at  this  location  is  determined  by  first  assuming  an  analytical  expression  for  the 
shock  shape.  8-1 


x  =  Rn  +  w  -  Rc  cot2p  [(1  + 


R2tan2p 

9 

R2 

c 


.1/2 


-  1] 


(B-9) 


Anderson,  J.  D.,  Hypersonic  and  High  Temperature  Gasdynamics,  McGraw-Hill  Book  Co.,  1989. 
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Rn  is  the  nose  radius  and  &)  is  the  shock  standoff  distance  given  by 

(o  3.24 

—  =0.143  exp  [  —  ]  (B-10) 

en  M* 

Rc  is  the  radius  of  curvature  of  the  shock  at  its  vertex  (assuming  a  hyperbolic 
shock  shape),  which  is  given  by 

=  1.143  exp  [  °'-4— I  (B-ll) 

rn  (M^-l)12 

P  is  the  limiting  shock  wave  angle  downstream  from  the  nose  curvature  region. 
By  differentiating  Equation  B-9  to  get  drs/dx  and  evaluating  it  at  rs,  the  shock  slope 
at  this  point  can  be  found.  Then,  the  change  in  entropy  across  the  shock  can  easily  be 
determined.  If  it  matches  the  assumed  value  at  the  body  point,  then  the  original 
assumption  was  correct  and  the  solution  can  proceed  to  the  next  body  station.  If  not, 
a  second  value  of  entropy  is  assumed  and  the  process  is  repeated.  If  *V  s  value  also 
fails  to  match,  then  a  variable  secant  iteration  process  can  be  started,  v»hi«,h  usually 
converges  quite  rapidly  to  the  correct  local  entropy. 

If  the  body  under  consideration  is  a  blunted  cone,  the  local  entropy  value  will  lie 
somewhere  between  that  found  behind  a  normal  shock  at  the  free-stream  Mach 
number  and  that  found  behind  the  attached  oblique  shock  associated  with  the 
solution  for  a  sharp  cone  of  the  same  half  angle.  Tables  B-l  and  B-2  show  normal 
shock  and  sharp-cone  entropy  values  for  several  cone  angles.  As  shown,  they  can  be 
substantially  different  with  the  difference  increasing  as  the  Mach  number  goes 
higher.  Tables  B-3,  B-4,  and  B-  5  illustrate  the  effects  of  using  the  variable  entropy 
approach  at  various  streamwise  stations  on  a  blunted  cone.  Values  are  given  for 
entropy,  pressure,  temperature,  density,  recovery  temperature,  and  heat  transfer 
rate  for  solutions  using  constant  sharp-cone  entropy,  constant  normal  shock  entropy, 
and  variable  entropy  for  the  windward  plane  of  a  15-deg  cone  at  Mach  10.6  and  angle 
of  attack  of  5  deg.  Pressures,  because  of  the  boundary  layer  assumption,  are  not 
affected,  but  as  can  be  seen,  temperatures  and  densities  can  vary  by  almost  a  factor  of 
three.  A  serendipitous  result  of  the  manner  in  which  the  primitive  variables  interact 
in  the  various  cases  is  the  relative  small  variation  of  heat  transfer  rates  for  the  three 
entropy  assumptions. 

Figure  B-2  gives  an  indication  of  how  the  entropy  varies  along  the  body  surface 
as  a  function  of  distance  downstream  from  the  stagnation  point  using  the  present 
approach.  As  can  be  seen,  at  a  distance  of  50  nose  radii  downstream,  the  entropy 
values  have  largely  stabilized  to  the  sharp-cone  value.  Figure  B-3  shows  the 
differences  in  heat  transfer  created  by  using  the  variable  entropy  approach  versus 
constant  normal  shock  entropy  for  the  case  of  the  windward  plane  on  a  15-deg  cone  at 
10-deg  angle  of  attack.  The  variable  entropy  solution  gives  slightly  higher  heat 
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transfer  values,  and  the  difference  increases  with  increasing  distance  from  the 
stagnation  point,  as  would  be  expected.  The  discrepancies  between  the  two  methods 
are  not  great,  typically  being  10  to  15  percent  at  most. 


TABLE  B-l.  NORMAL  SHOCK  AND 
CONE  ENTROPIES,  M  =  8, 
ALTITUDE  =  150,000  FT 


Cone  Half  Angle 
(deg) 

Normal  Shock 
Entropy 
(ft2/sec2-°R) 

Cone  Entropy 
(ft2/sec2-°R) 

5 

60631 

51914 

10 

60631 

52206 

15 

60631 

52924 

20 

60631 

53917 

25 

60631 

54963 

30 

60631 

55962  ! 

TABLE  B-2.  NORMAL  SHOCK  AND 
CONE  ENTROPIES,  M  =  15, 
ALTITUDE  =  150,000  FT 


Cone  Half  Angle 
(deg) 

Normal  Shock 
Entropy 
(ft2/sec2-°R) 

Cone  Entropy 
(ft2/sec2-°R) 

5 

69548 

52139 

10 

69548 

53706 

15 

69548 

55766 

20 

69548 

57632 

25 

69548 

59374 

30 

69548 

61391 
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TABLE  B-3.  ENTROPY  EFFECTS  ON  15-DEG  CONE,  X  =  0.528  FT, 
M  =  10.6,  ALTITUDE  =  150,000  FT,  ANGLE  OF  ATTACK  =  ODEG, 
WINDWARD  PLANE 


Parameter 

Cone  Entropy 

Stagnation 

Entropy 

Variable 

Entropy 

Entropy  (ft2/sec2-°R) 

45394 

52308 

51444 

Pressure  (lb/ft2) 

53.78 

53.74 

53.74 

Temperature  (°R) 

376 

1190 

1036 

Density  (slug/ft3) 

0.000083229 

0.000025946 

0.000029760 

Recovery  temperature  (°R) 

1763 

1854 

1830 

Heat  transfer  (BTU/ft2-sec) 

3.0 

2.5 

2.6 

TABLE  B-4.  ENTROPY  EFFECTS  ON  15-DEG  CONE,  X  =  1.025  FT, 
M  =  10.6,  ALTITUDE  =  150,000  FT,  ANGLE  OF  ATTACK  =  ODEG, 
WINDWARD  PLANE 


Parameter 

Cone  Entropy 

Stagnation 

Entropy 

Variable 

Entropy 

Entropy  (ft2/sec2-°R) 

45394 

52308 

50605 

Pressure  (lb/ft2) 

53.78 

53.74 

53.74 

Temperature  (°R) 

376 

1190 

906 

Density  (slug/ft3) 

0.000083229 

0.000025946 

0.000034120 

Recovery  temperature  (°R) 

1763 

1854 

1806 

Heat  transfer  (BTU/ft2-sec) 

2.1 

1.8 

1.9 

TABLE  B-5.  ENTROPY  EFFECTS  ON  15-DEG  CONE,  X  =  2.033  FT, 
M  =  10.6,  ALTITUDE  =  150,000  FT,  ANGLE  OF  ATTACK  =  0  DEG, 
WINDWARD  PLANE 


Parameter 

Cone  Entropy 

Stagnation 

Entropy 

Variable 

Entropy 

Entropy  (ft2/sec2-°R) 

45394 

52308 

49135 

Pressure  (lb/ft2) 

53.78 

53.74 

53.74 

Temperature  (°R) 

376 

1190 

429 

Density  (slug/ft3) 

0.000083229 

0.000025946 

0.000072911 

Recovery  temperature  (°R) 

1831 

1854 

1829 

Heat  transfer  (BTU/ft2-sec) 

1.4 

1.3 

1.4 
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FIGURE  B-2.  VARIATION  OF  ENTROPY  ALONG  WINDWARD  SYMMETRY  PLANE 
OF  1.5-IN.  NOSE  RADIUS,  15-DEG  HALF-ANGLE  CONE 


X/Rn 


FIGURE  B-3.  EFFECT  OF  ENTROPY  ON  HEATING  RATES  ON  WINDWARD 
SYMMETRY  PLANE  OF  1.1-IN.  NOSE  RADIUS,  15-DEG  HALF-ANGLE  CONE 
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